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Through 2D granular Couette flow experiments, we probe failure and deformation of disordered 
solids under shear. Shear produces smooth affine deformations in such a solid and also irresversible 
so-called non-affme particle displacements. We examine both processes. We show that the non- 
affine part is associated with diffusion, and also can be used to define a granular temperature. 
Distributions for single particle non-affme displacements, Sn, satisfy Pi(Sn) oc exp[— |<5ri/Ar| a ] 
(q ~ 2). We suggest that the shear band forms due to a radially outward diffusive flux/non- affine 
motion which is balanced in the steady state by inward diffusion due to density gradients. 



The way in which disordered solids undergo plastic i.e. 
irreversible deformation is a topic of considerable recent 
interest 0, S B 0, IS @, 0, H, in part, because the phe- 
nomena appear to be general across a large collection of 
materials that include glasses, foams, colloids and gran- 
ular materials, among others. Here, we focus on granu- 
lar materials. Despite a variety of experiments on bulk 
deformation of granular materials, the microscopic ba- 
sis for plasticity remains elusive. Molecular dynamics 
simulations @, Si,!,!] have provided some insight, al- 
though a number of these studies involve particles that 
interact via frictionless forces. At this point, there is rel- 
atively little experimental work that directly addresses 
the connection between the microscopic and larger scale 
nature of granular plasticity. The goal of the present 
work is to begin to fill this void. 

There is, however, a relatively extensive literature that 
considers microscopic plastic deformation in amorphous 
molecular solids [1, H, In the limit of low temper- 

atures where thermal noise is irrelevant, it is expected 
that the plastic behavior in such materials is similar to 
what occurs in a granular materials. Both microscopic 
and granular materials exist in disordered jammed states, 
and when sheared, undergo irreversible configurational 
changes. In the granular case, the configurational struc- 
tures also involve force chains, long filimentary structures 
that form so as to resist shear, e.g. Fig. [TJ 

In recent theoretical work, Falk and Langer (FL) 1] 
developed models of elasto-plastic deformation of amor- 
phous solids which have also been applied recently to 
granular materials 0, Q • FL consider localized regions 
that are susceptible to irreversible failure. Associated 
with any smooth affine structural change or deformation 
in such a region, there is a non-affine part which leads 
to irreversibility/plasticity. In their analysis, FL [l| con- 
sider small clusters of N particles. They fit local deforma- 
tions occurring over a time At to an affine deformation. 
To identify localized plastic rearrangements, they define 
a quantity (D^ m ) which characterizes departures from 
the local affine deformation. 

In this letter, we characterize the motion generated in 
2D granular Couette (shear) flow in terms of affine and 
non-affine components. For this system, deformation is 




FIG. I: Time sequence of images emphasizing the spatially 
inhomogeneous and temporally intermittent nature of force 
chains (bright lines) and active areas of non-affine deforma- 
tion (shaded disks). Each disk shows a local instantaneous 
value for D^ in (see text), where brighter corresponds to larger 
Dm iri . Images are separated in time by 5 At (see text), and 
run left-to-right, top-to-bottom. The shearing wheel, located 
near the bottom of each image, rotates counterclockwise. Nor- 
mally black, the edge of the wheel is highlighted for visability. 



localized to a shear band of width less than 10 grain diam- 
eters. Despite the relatively rapid variation of quantities 
with radial distance, it is still possible to characterize the 
local affine and non-affine components of the flow. How- 
ever, there are some differences, discussed below, between 
the case of homogeneous pure shear studied in FL, and 
the present case, which is spatially varying and also in- 
volves rotation of small volume elements. 

The experimental apparatus used here and described 
in greater detail elsewhere 0, is a 2D Couette shearing 
experiment. A dense packing of ~ 45, 000 bidisperse 
photoelastic disks lie on a flat horizontal surface. In 
that plane, the disks are bounded by an inner shear- 
ing wheel and an outer concentric ring. The wheel ro- 
tates slowly (/ = inverse rotation period = lmHz in the 
present experiments) to impose quasistatic shear. Parti- 
cles are marked with lines for easier tracking of position 
and rotation. The apparatus is illuminated from below 
and imaged from above by two cameras. One camera 
observes the grains through crossed polarizers while the 
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other simultaneously takes unpolarized images. In this 
way, we record both stresses (through photoelastic mea- 
surement with polarized images) and kinematic informa- 
tion (by particle tracking with unpolarized images). We 
note that the motion is confined to a shear band near 
the shearing wheel. The tangential velocity v(r) and lo- 
cal shear rate dv/dr decrease faster than exponentially 
with radial distance from the shearing surface and have 
fallen by three orders of magnitude from their value at the 
shearing wheel for r ~ 8<i[l(|. Here, we choose r = to 
correspond to the surface of the shearing wheel. In addi- 
tion to the radial velocity gradient, there is also a density 
gradient, i.e. Reynolds dilatancy. The density /packing 
fraction near the wheel can be lower by as much as 2/3 
of the density far away. 

In order to characterize the motion, we consider lo- 
calized regions falling within 2.2<i of a reference parti- 
cle, where d is the mean particle diameter. (FL use a 
similar sampling radius of 2.5asL where asL is an inter- 
atomic distance.) Typically, we find N ~ 15 particles 
in such a neighborhood, including the particle of inter- 
est. We then follow the motion of these N particles over 
a time At = 1.85s that is short compared to l/f. We 
are interested in the deformation of this patch. To this 
end, we determine the location of each particle relative 
to the patch center of mass (CoM), both initially (r^), 
and after the small elapsed time (r^). We determine the 
smooth mapping that best describes the affine deforma- 
tion, by least squares fitting data for the r t and r\ in a 
cluster to the form r\ = Er,, where E is a 2 x 2 ma- 
trix. D'^ lin — E(r- — Erj) 2 is the parameter of Falk and 
Langer, with the slight modification that we must con- 
sider motion relative to the CoM, since there is net flow 
in our system. E is in general not symmetric, and hence 
its eigenvalues may be (and often are) complex. This is 
a consequence of the typical deformation in the present 
flows: on average, there is both deformation of a volume 
element, characterized by a symmetric tensor, F, and 
rotation, characterized by a rotation tensor, R#. This 
differs from the case of pure shear, where there is no ro- 
tation on average. We write E = RqF without serious 
ambiguity. For instance, inverting the order of F and 
R# does not change 9 or the eigenvalues of F. Also, 9 
is uniquely determined if we restrict — tt/2 < 9 < n/2, 
which is not an issue for the small deformations consid- 
ered here. We can then express the deformation in terms 
of the eigenvalues, e,, of e = F — I, where e is the conven- 
tional strain tensor, and I is the unit tensor. Together 
with 6, they characterize the affine part of the local evo- 
lution. D^m characterizes the extent of the nonaffinc 
motion for the small region in question. We note that all 
measures of the deformation, both affine and non-affinc 
are spatially and temporally intermittent, as seen for in- 
stance in Fig. [TJ Hence, it is important to consider the 
probability distribution functions (PDF's) of these mea- 
sures as functions of r. 

We first consider data for the affine deformations. In 
Fig. [5] we show PDF's for the eigenvalues of e, which we 
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Affine Deformation, 8e = £., - ^ l 

FIG. 2: PDF's of the affine dilation, ei + £2 (top), and affine 
shear deformation, 62 — ei (bottom). Each series corresponds 
to a radial bin of width d. The notation r = 1 — 2d, etc. 
indicates that the central particle of the cluster was located 
radially between d and 2d from the wheel, etc. 



give in the form 5e = 62 — e± and 2e = ei+€2- The different 
data sets in each figure are binned over radial widths 
of size d, and correspond to a series of radial distances, 
each d apart, from the shearing wheel. The first series is 
centered at d/2 from the shearing wheel. We expect that 
e has zero mean in the steady state, since it corresponds 
to the dilation of the local patches; 5e has non-zero mean 
and a characteristic shape that reappears in the data for 
£)^ in , discussed below. PDF's for 9 resemble those for 
<5e; mean values for 9 are typically a few degrees or less. 

We plot the PDF's of -D^ in binned at different dis- 
tances r from the shearing surface in Fig. OUtop). The 
peak in the PDF of -D^ in tends to be at larger values of 
Dmin f° r r closer to the shearing wheel, although the po- 
sition of the peak does not vary monotonically with r, a 
point that we consider below. The individual PDF's are 
well fitted by the form f(x) — Cix C2 e~ x ^ C3 , which we 
discuss further below. However, the exponent C2 and the 
scaling factor C3 in the exponential vary from fit to fit. 
In Fig. [3{bottom) , each of the data sets is rescaled in the 
vertical direction by the peak fit value and in the horizon- 
tal direction by the mean (D^ in ). We see that the data 
at different r collapse fairly well onto a single curve. The 
averages (-D^ in ) (inset), are set by the local shear rate 
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FIG. 3: Top: PDF's of D 2 min (units of d 2 ) at various distances 
r from the shearing surface (in units of d) for 1080 images. 
Bottom: Rescaled data from top. Each data set is fitted and 
rescaled by the peak magnitude and mean value of fl^ jn . 
Inset: (fl„ in ) vs. local shear rate. Note the drop in {D^„) 
(departure from red line) at high 7, i.e. near wheel. 



within the shear band (inset, bottom). We previously 
observed a similar dependence for diffusivity, D, on local 
shear rate 7 , i.e. D oc 7 Q. In Fig. 01 we plot {D^ in ) and 
radial and tangential diffusivities (D rr and Dee) rescaled 
such that all quantities are 1 at r = 2d. The diffusivity is 
therefore proportional to the mean measure of local plas- 
tic rearrangements (D^ in ). It is perhaps not surprising 
that diffusivity and D^ in are related, as they both mea- 
sure non-afhne displacements. Dimensionally, D 2 nin has 
units of x 2 and is measured over a specified time interval 
At. Diffusivity D is set by D ~ ((Ax) 2 } /At. In these ex- 
periments, displacements scale by some fraction of d over 

e.g. 



times scales of order the local inverse shear rate 7 
) 2 ■ 

nit n 



D ~ 7 d 2 and D 2 mm ~ 7 Atd 2 . 

The decrease in diffusivity and (D^J for r w in 
Fig. |4] results from several effects. Particles in contact 
with the shearing surface are dragged quasi-ballistically 
in the azimuthal direction Q. In addition, radial diffusion 
is limited to the outward direction next to the shearing 
wheel. Roughly, one might expect a decrease by 1/2 for 
particle variances for radial displacements (and hence ra- 
dial diffusivity, D rr ) near the wheel. 

D 2 nin is a mesosopic property, i.e. a locally averaged 
measure of a grain-scale process. The scaling property 
observed for the PDF's of D J^ in suggest that there may 
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FIG. 4: Rescaled (D^ in ), widths of non-affine PDF's, Pi(Sri), 
diffusivities, and peak values of P(9) and P(e2 — ei) versus 
radial distance from the shearing wheel. All quantities are 
scaled to 1 at r = 2d. For reference, at r — 2d, 7 = 0.60s -1 , 
D rr = 0.0054d 2 /s, Deo = 0.011d 2 /s, < D 2 min >= 0.38d 2 , 
Opeak = 0.045rad, (e2 — £i) P eak = 0.047d. Square widths for 

other PDF's are: P(e 2 - ei)-0.0016d 2 , (Pi(5n) radial) 2 - 

0.0050d 2 , (Pi(Sn) - - azimuthal)-0. 0053d 2 . 



be an underlying microscopic phenomena of interest. To 
this end, we examine in Fig [5] (top), the PDF's Pi(5r\) 
of the individual non-affine particle displacements, <5r.; = 
r\ — Erj. The PDF for P\ is approximately gaussian, 
Pi(Sri) ~ Aexp[-(<5r l /Ar) Q ], where a < 2. To see this, 
we plot in Fig. \E\ (middle), log\log(Pi(8r i ))\ as a func- 
tion of log(Sri), with Pi normalized so that P\(0) = 1. 
These data are striking in as much as they resemble the 
velocity distributions for granular gas-like states pdl. [l~2l|. 
even though the basic kinetic theory assumption of short- 
lived uncorrelated collisions/contacts is manifestly vio- 
lated. Interestingly, PDF's for the total displacements, 
r\ — ri, are exponentially distributed, Fig [5] (bottom) . 
D 2 min is constructed by choosing N particles in a clus- 



ter, and then computing D 



2 

min 



{Z N 6r 2 ). The PDF for 



Df nin , constructed from clusters of N particles is 



P(DLn)= / PN(6r 1 ,..5r N )6(D 2 mm -Z(6r l ) 2 )d6r 1 ...d5r N , 



(1) 

where Pjv gives the PDF for the N non-affine particle 
displacements. We do not a priori know Pn(5vi, ..5rjv)- 
However, if we assume that the individual particle 
displacements are uncorrelated and given by gaussian 
PDF's: Pi(Sn) = Bexp[-(Sri/Ar) 2 }, then the above 
integral for the PDF for D 2 in follows with a bit of al- 
gebra. (Although neither of these assumptions is rig- 
orously true, neither is markedly unreasonable.) Then, 
P(D 2 min ) = C 1 (^ 2 am ) JV - 1 > exp(-^ m /C 2 ), where the a 
are constants. This form is qualitatively consistent with 
the experimentally obtained PDF's for D 2 min (Fig. [3]). 
Applying this strictly is probably not warranted, since 
some correlation is likely, and the measured P\ 's are not 
truly gaussians. 

The widths of the P\(5ri) (Fig. can be taken as an 
effective granular temperature associated with the non- 
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FIG. 5: Data for individual particle non-affine radial displace- 
ments, Pi(Sri). Top: on linear scales, Middle: on loglog-log 
scales (solid black line shows slope 3/2). PDF's for the non- 
affine azimuthal displacements are similar. Bottom: PDF's 
of total displacements, r[ — n, on log-lin scales. 



affine/diffusive displacements of the particles (hence pre- 
sumably a configurational measure). However, we give a 
note of caution here concerning the relation of this tem- 
perature to the overall energy in the system. The elastic 
energy and elastic energy fluctuations for this system are 
several orders of magnitude larger than the correspond- 
ing kinetic energies associated with both the affine and 
non-affine motion. In addition, the motion considered 
here is driven by the wheel displacement, not its veloc- 
ity. If the experiments were slowed by any amount, the 
basic phenomena would appear essentially the same when 
expressed in terms of displacements. 

In closing we offer the following insight into the for- 
mation of a shear band from an initially homogeneous 
state. Both the affine and non-affine motion are driven 
by and extract energy from the wheel displacement and 
the stored elastic energy. Non-affine deformation is asso- 
ciated with diffusive motion, whose source is the stirring 
from the wheel. When the wheel begins to turn, there 
is a radially outward flux driven by the excitation of the 
shearing wheel (there can be no inward diffusion there). 
This flux is the source of the (Reynolds) dilation next to 
the wheel. Such a flux is diffusive in the steady state, 
and is then counterbalanced by the resulting gradients 
of density/free volume which cause a compensating re- 
verse mass flux. This state is also mechanically stable 
due to the reduction in density and contact number near 
the wheel. The system is less resistant to shear near the 
wheel due to contact number depletion, leading to local- 
ization of motion. In the present system, which has a 
fixed volume, one expects an increase in density outside 
the dilated region, which tends to limit the growth of 
the shear band, in contrast to recent 3D experiments [P3T] 
with an open top where the growth of the shear band 
is much less constricted. An interesting question is the 
extent to which similar phenomena occur in other disor- 
dered solids. 
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